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Abstract 

We study via the renormalization group equations at one-loop order the perturbativity and 
vacuum stability of a two-singlet model of cold dark matter (DM) that consists in extending the 



(— I ' Standard Model with two real gauge-singlet scalar fields. We then investigate the regions in the 

parameter space in which the model is viable. For this, we require the model to reproduce the 



observed DM relic density abundance, to comply with the measured XENON 100 direct-detection 
upper bounds, and to be consistent with the RGE perturbativity and vacuum-stability criteria up 



^ . to 40TeV. For small mixing angle 9 between the physical Higgs h and auxiliary field, and DM-/i 

C^ ■ 

mutual coupling constant Aq , we find that the auxiliary-field mass is confined to the interval 

^^ ' 116GeV — 138GeV while the DM mass is mainly confined to the region above 57GeV. Increasing 

"^ ' (4) 

^D ■ enriches the existing viability regions without relocating them, while increasing X\ shrinks 

m 

them with a tiny relocation. We show that the model is consistent with the recent Higgs boson-like 
-^ discovery by the ATLAS and CMS experiments, while very light dark matter (masses below 5GeV) 

X' 

H ' is ruled out by the same experiments. 
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I. INTRODUCTION 

Now that there is more and more compelhng evidence that the discovery in the ATLAS 
and CMS experiments at the LHC is a Higgs particle with a mass rrih — 125GeV l|, |2|, 
one of the main focuses of particle physics is an understanding of the still elusive nature of 
dark matter (DM), believed to account for about 26% of the energy content of the universe 
|3|. Side by side with observation, models beyond the Standard Model (SM) are devised 
to account for weakly interacting massive particles (WIMPs) as plausible candidates for 
dark matter. Such models range from the sophisticated ones, those bearing intricate under- 
lying symmetries and mechanisms of symmetry breaking, to the more simpler ones, those 
extending the SM without any particular assumption regarding a deeper structure. 

Since the findings at the LHC have so far yielded no particular clue as to possible particle 
structures beyond those of the Standard Model, it is therefore still consistent to try to 
recognize dark matter as simple WIMPs extending the Standard Model, with no further 
assumptions as to inner structures. The simplest of such extensions is a one real electroweak- 
singlet scalar field, interacting with visible matter only via the SM Higgs particle, first 
proposed by Silveira and Zee [^ and further studied in |5|-l7|. In this minimal model, with 
DM masses lighter than lOGeV, the Higgs would be mainly invisible, which is excluded by 
the recent measurement of the Higgs signal at the LHC [8]. Also, for DM masses in the 
lange 7GeV - 60GeV, the model is ruled out by the data from XENON 10 [9] and CDMS 



II IDI], except around the resonance 62GeV at which the Higgs-DM coupling is extremely 



small. Furthermore, DM masses between 65GeV and 80GeV are excluded in this one-singlet 



extension by the XENON 100 experiment 11 



a complex scalar singlet extension of the SM 12 1 



Note that similar conclusions hold also for 



So, given the difficulties this minimal model has with existing experimental and observa- 
tional DM data p^] , we have proposed in [ij] an extension to the Standard Model with two 
real electroweak singlets, one being stable, the DM candidate, and the other an auxiliary 
field with a spontaneously broken Z2 symmetry. Based on the DM relic-density and WIMP 
direct detection studies, we have concluded that this two-singlet model is capable of bearing 
dark matter in a large region of the parameter space. Further constraints on the model 
as well as some of its phenomenological implications have been studied in p^] where rare 
meson decays and Higgs production channels have been considered. 



In the present work, we further the study of the two-singlet model and ask how high in 
the energy scale it is computationally reliable. A standard treatment is the investigation of 
the running of the coupling constants in terms of the mass scale A via the renorrmalization 
group equations (RGE). We believe one-loop calculations are amply sufficient for the present 
task; higher loops could be considered if the situation changes. 

The two standard issues to monitor are the perturbativity of the scalar coupling constants 



and the vacuum stability of the theory. These issues were studied in 16| for the complex 
scalar singlet extension of the SM, and it was shown that the vacuum-stability requirement 
can affect the DM relic density. Specific results from such studies depend on the cutoff scale 
Am of the theory. Reversely, imposing perturbativity and vacuum stability may indicate at 
what Krri the two-singlet extension is valid. In the early parts of this work, the second point 
of view is adopted, whereas in the later part, the ffist is taken. Furthermore, up to only 
recently, it has been anticipated that new physics such as supersymmetry would appear at 
the LHC at the scale A ~ ITeV. Present results from ATLAS and CMS indicate no such 
signs yet. One consequence of this is that the cutoff scale A^ naay be higher. As we shall 
discuss, we find that it can be ~ 40TeV. 

This work is organized as follows. After this introduction, we recapitulate in section II 
the essentials of the two-singlet model necessary for the RGE calculations. In section III, we 
discuss the running of the scalar coupling constants when we switch off the non-Higgs SM 
particles. This gives us a ffist understanding of how high in the energy scale perturbativity 
is allowed. It helps also seeing the subsequent effects of the other SM particles. Section IV 
discusses the full RGEs. Vacuum stability gets into focus with the Higgs coupling constant 
turning negative at some scale. The mass scales at which non-perturbativity and non- 
stability set in are different, and so a choice for A^ has to be made. Section V attempts at 
finding the regions in the parameter space in which the model is predictive. In addition to 
the DM relic-density constraint and the perturbativity-stability criteria deduced from the 
previous two sections, we impose on the model to be within the current direct-detection 
experimental bounds. Section VI is devoted to concluding remarks. 



II. THE TWO-SINGLET MODEL 

The model is obtained by adding to the Standard Model two real, spinless, and Z2- 
symmetric SM-gauge- singlet fields. One is the dark matter field Sq with unbroken Z2 sym- 
metry, and the other an auxiliary field Xi with spontaneously broken Z2 symmetry. Both 
fields interact with the SM particles via the Higgs doublet H. Using the same notation as 



in IJ], the potential function that involves Sq, H and Xi is: 



+|S„^ + ^ (WHf + |xt + ^SlH^H + 'JfSlxl + ^H^Hxl (2.1) 

where mg, /i^and ^\ and all the coupling constants are real positive numbers^. 

We are interested in monitoring the running of the scalar coupling constants. A one-loop 
renormalization-group calculation yields the following /3-functions: 

/5a = ^ (^A^ + \l + \l- 48A^ + 8AA? - 3A/ - \g'' + ^g^g'^ + \g' 

/5r,oi = Y^ (4^01 + ^o??oi + TliTlm + 4AoAi) ; 

/3ao = 77^ ( 4Ag + Aor/o + 2AoA + r^oiAi + 12AoA2 - ^X^g^ - ^X^g'^ 



1 / 9 3 

^— ( 4A? + Air^i + 2AiA + r/oiAo + 12AiA? - -Ai/ - -, 



/5ai = TTT^ ( 4A? + Air^i + 2AiA + r/oiAo + 12AiA? - -Ai^?^ - -Ai^?'^ ) . (2.2) 



As usual, Pg = dg/dlnA where A is the running mass scale, starting from Aq = lOOGeV. 
The constants g, g' and gg are the SM and strong gauge couplings, known [l7| and given to 
one-loop order by the expression: 

G (A) = ^^^°^ , (2.3) 



1 - 2aG G^ (Ao) In ( A) 



where ac = gg^, ^, j^, and G (Aq) = 0.65, 0.36, 1.2 for G = g,g',gs respectively. The 
coupling constant At is that between the Higgs field and the top quark. To one-loop order, 



^ The mutual couplings can be negative as discussed below, see (13.1 



it runs according to |17] : 

with Xt (Ao) = "'^^ = 0.7, where v is the Higgs vacuum expectation value (vev) and 
rrit the top mass. Note that we are taking into consideration the fact that the top-quark 
contribution is dominant over that of the other fermions of the Standard Model. 

The model undergoes two spontaneous breakings of symmetry: one of the electroweak, 
with a vev v = 246 GeV, and one of the Z2 symmetry (xi field), with a vev Vi we take in 
this work equal to 150 GeV. Above v, the fields and parameters of the theory are those 
of (12. ip . Below vi, the (scalar) physical fields are 5*0 (DM), h (Higgs) and S-\ (auxiliary), 



with parameters (masses and coupling constants) given in Eqs. (2.2-2.15) of IJ]. We take 
the values of the physical parameters at the mass scale Aq = 100 GeV. There are originally 
nine free physical parameters. The two vevs v and vi are fixed, as well as the mass of the 



physical Higgs field nih = 125 GeV 



(4) 



21]. Also, the physical mutual coupling constant rj^^ 



between Sq and Si is determined by the DM relic-density constraint 18|, which translates 
into the condition: 

(^i2^ann)- 1.7 X 10-9 GeV-^ (2.5) 

where (fi20"ann) is the thermally averaged annihilation cross-section of a pair of two DM 
particles times their relative speed in the center-of-mass reference frame. This constraint is 
imposed throughout this work, together with the perturbativity restriction < r/g^ < v^ivr 
on its solution. The remaining free parameters of the model are the physical mutual coupling 
constant Aq between h and 5*0, the mixing angle 9 between h and 5*1, the DM mass uiq, 
the mass mi of the auxiliary physical field Si, and the DM self-coupl ing constant t]q. This 



fl 



latter has so far been decoupled from the other coupling constants [Ij, Il5| , but not anymore 
in view of (12. 2p now that the running is the focus. However, its initial value rjQ (Aq) is 
arbitrary and its /3-function is always positive. This means r]o (A) will only increase as 
A increases, quickly if starting from a rather large initial value, slowly if not. Therefore, 
without loosing generality in the subsequent discussion, we fix rjQ (Aq) = 1. Hence, here too 
we still effectively have four free parameters: Ai , 6, mo, and mi. The initial conditions for 



the coupling constants in (12. ip in terms of these physical free parameters are as follows: 



Vi (Ao) 

A(Ao) 

Ai (Ao) 

'701 (Aq) 

Ao (Ao) 
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nit. 
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Till + ml — Iml — m^l (cos (29) 
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sin {29) 



r/ii 



m 



h\ J 



r^^^ cosH - XC^ sm^ 9 



(2.6) 



cos {29) 

Note that normally, as we go down the mass scale, we should seam quantities in steps: at 
f, fi, and Aq. However, the corrections to (12.61) are of one-loop order times In ;^ or In^, 
small enough for our present purposes to neglect. 



III. SCALARS ONLY 

To see the effects of the scalar couplings only and how up in the mass scale the model 
can go, we switch off the non-Higgs SM couplings in (12. 2p . The perturbativity constraint 
we impose on all dimensionless scalar coupling constants is G (A) < yAnx. Vacuum stability 
means that G (A) > for the self-coupling constants r^o. A, and r^i, and the conditions: 



1 ^— , r— 1 , , — 

-yriaX < Ao < V47r; --VWi < ^701 < v47r; 



-7:\A?iA < Ai < 
6 



47r (3.1) 



for the mutual couplings Ao, "^oi, and Ai. Also, as a start, we let the masses mo and mi vary 
in the interval IGeV - 200GeV. 

Fig. [1] displays a typical running of the scalar self-coupling constants, from Aq = lO^GeV 
up to lO^^GeV. As is expected for scalars only, all coupling constants are increasing functions 
of the scale A, with different but increasing rates. Also, the larger value the coupling starts 
from at Aq, the faster it will go up. Fig. |2] shows the running of the mutual coupling constants 
for the same values of the parameters. For these values, the mutual coupling constants start 
well below 1, and so run low; they are very much dominated by the self- couplings. This 
situation will stay for Aq and Ai in all regions, but not for r^oi. 

The first coupling constant that leaves the perturbativity bound \/Atx is r^i, the self- 
coupling constant of the auxiliary scalar field Xii ^^ about 1260 TeV for this set of values of 
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Ao*"^' = 0.01, e = 10°, OTo = 55GeV, mj = 1 lOGeV 
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FIG. 1: The running of the self-couphngs (scalars only). The self-couphng r]i of the auxihary field 
Xi dominates over the Higgs self-coupling A. 

the parameters. This behavior is in fact typical. Indeed, rji starts above 2 at Aq in all the 
parameter space, much higher than all the other coupling constants - only rjoi can compete 
with it in some regions. As it intervenes squared in its own /3-function, it will also move 
up quicker. More precisely, from (12. 6p . we see that ?7i (Aq) depends on mi and 9 only. The 
effect of the mixing angle 6 is small. As a function of mi, starting from about 2, r/i (Aq) 
decreases slightly until rrih and then picks up. It will pass the perturbativity bound V^ at 
about mi ~ 160GeV. This means that the region mi > 160GeV is automatically excluded 
from the outset. In actual situations, given the positive-slope RG running of rji (A) and even 
in the case of the full RGE (see below), perturbativity puts a stricter upper bound on mi, 
irrespective of the other parameters of the model. 

In Fig. [H the Higgs self-coupling A starts just above 0.6 and does not pick up much when 
running. This behavior is typical too. Indeed, A (Aq) is also a function of m-i and 6 only, 
see (12. 6p . For a given 6, it will increase as a function of mi to reach 3 {rrih/v) at mi = rrih, 
equal here to 0.77 for ruh = 125GeV. Then it continues to increase, but with a smaller 
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FIG. 2: Running of the mutual couplings (scalars only). For these values of the parameters, all 
three are well below the self-couplings. 

slope. The mixing angle 6 enhances the behavior of A (Aq) as a function of mi, but for, say^ 
6 = 15°, A (Aq) will be less than 0.85 at rrii = 160 GeV. This situation implies that when 
running as a function of the scale A, the Higgs self-coupling A will increase, but will hardly 
reach 1 before, say r]i, leaves the perturbativity bound. 

Increasing the dark-matter mass does affect the running of the couplings. Figs. [3] (self- 
couplings) and m (mutual couplings) display such effects for ttiq = lOOGeV. Among the 
self-couphngs, rji is still dominant, but tailed more closely by rjQ this time. For both, the 
positive acceleration is accentuated, something that makes rji leave the perturbativity region 
much earlier, at about 6.3TeV. By contrast, the running of the Higgs self-coupling A stays 
flat. 

The major effect of increasing tuq is on the mutual coupling rjoi, between the DM field 5*0 
and the auxiliary field Xi- Indeed, in Fig. [2] where rriQ = 55GeV, r/oi started and ran small 



^ We are implicitly confining the mixing angle to small values, a situation inferred from our work 15| on 
the phenomenological implications of the model. This is discussed later in section V. 
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Ao*'" = 0.01, e = 10°, OTo = lOOGeV, otj = llOGeV 
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FIG. 3: Running of the self-couplings (scalars only) for a larger DM mass, rji is still dominant, 
even if r/o picks up faster behind it. 

like the other two mutual couplings. Here, whereas Aq and Ai (both Higgs related) stay close 
to zero, rjoi starts above 2.5 and runs up fast. In fact, for these values of the parameters, it 
leaves the perturbativity region earlier than rji, at about 2.5TeV. 

Increasing the auxiliary- field mass rrii has a similar effect: it enhances the positive ac- 
celeration of the self-couplings rji and tjq while leaving A flat, and boosts up the mutual 
coupling ?7oi away from Aq and Ai, which both remain not far from zero. It also makes t]i 
and ?7oi leave the perturbativity region earlier, without rjoi necessarily taking over from t]i. 

Increasing Aq has also an effect. Figs. (self) and El (mutual) show the running for 



A 



(4) 



0.4. The self-coupling rji dominates and leaves the perturbativity region at about 



251TeV. The mutual coupling r/oi is raised above 1 at Aq and so can run high, while Aq and 
Ai stay here too just above zero. It leaves the perturbativity region at about 794GeV, well 
behind %. Higher values of Ag are more difficult to achieve as the relic-density constraint 
(12. 5p may not be satisfied |l4l |. 

Finally, changing the mixing angle 6 has little effect on the self-coupling constants. It 
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FIG. 4: Running of the mutual couplings (scalars only) for a larger DM mass, ryoi starts above 1, 
much higher than the two others, even higher than the self-coupling ri\. 

helps the mutual coupling constants r^oi and Ai start higher, but not by much: they stay 
with Aq well below one. 

IV. THE FULL RGE 

In the previous situation, 'scalars only', all running coupling constants were positive, 
and so there were no issues related to vacuum stability. We now reintroduce the other SM 
particles and see their effects. Fig. [7] displays the behavior of the self-couplings under the 
full RGE for the same values of the parameters as in Fig. [1] (scalars only). The dramatic 
effect is on the Higgs self-coupling constant A which quickly gets into negative territory, at 
about 15TeV, thus rendering the theory unstable beyond this mass scale. This is better 
displayed in Fig. |H] where the RG behavior of A is shown by itself. Such a negative slope for 
A is expected, given the negative contributions to fi\ in fl2.2p . Here too ?7i is dominant over 
the other couplings and still controls perturbativity, leaving its region much later, at about 
1600 TeV, farther away from the situation 'scalars only'. This looks to be a somewhat general 
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FIG. 5: Running of the self-couplings (scalars only) with a larger physical Higgs self-coupling Aq . 
The self-coupling r]\ is still dominant. 

trend: the non-Higgs SM particles seem to flatten the runnings of the scalar couplings. 

The runnings of the mutual coupling constants for the same set of parameters' values is 
displayed in Fig. O They also get flattened by the other SM particles, but they stay positive. 
Here too they dwell well below the self-couplings, as in the 'scalars only' case. In fact, many 
of the effects on the running coupling constants coming from varying the parameters are 
similar to those of the previous situation since the SM particles do not intervene in the 
initial values of the couplings (self and mutual) at Aq. This means that increasing vtiq and 
vtix will raise the mutual coupling ?7oi and not the two others, higher than ?7i in some regions. 
For example. Fig. [TUJ shows the running of the self-couplings when mo = lOOGeV. Both r/i 
and r/o run faster but A is little affected. Fig. [11] shows the running of the mutual couplings 
from the full RGE also at tuq = lOOGeV. As in the case 'scalars only', larger ruQ boosts up 
riQi (Aq), much higher than Aq and Ai, at about 2.2 here, which makes it run quickly high, 
leaving the perturbativity region before rji, as in the case 'scalars only'. 

Raising Aq will also make the self-couplings ?7i and r]Q run faster while affecting very 
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FIG. 6: Running of mutual couplings (scalars only). Larger Ag helps r/oi rise well above Aq and 
Ai, but not enough to win over the self-coupling r]i. 

little A. It will also make the mutual coupling rjoi starts higher, and so demarked from Aq 
and Ai. By contrast, the effect of 6 is not very dramatic: the self-couplings are not much 
affected and the mutuals only evolve differently, without any particular boosting of rjoi. 

V. REGIONS OF VIABILITY 

The foregoing discussion showed us how the scalar parameters of the two-singlet model 
behave as a function of the mass scale A. From the situation 'scalars only' we understood 
that the two couplings that control perturbativity are t]i and t]qi. The full RGE brought in 
stability: the change of sign of A is the vacuum stability criterion to use. Equipped with 
these indicators, we can try to investigate in a more systematic way the viability regions of 
the model, regions in the space of parameters in which the model is predictive. Remember 
that this model has four parameters: the dark-matter mass rriQ, the physical auxiliary field 
mass mi, the physical Higgs self-coupling Aq , and the mixing angle 9 between the physical 
Higgs and the auxiliary field. The way we proceed is to vary Aq and 6 and try to find the 
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FIG. 7: Running of the self-couplings (full RGE). ryi controls perturbativity and the Higgs coupling 
A becomes negative quickly. 

regions of viability of the model in the (mo, mi)-plane. 

We have by now a number of tools at our disposal. First the DM relic-density constraint 
(I2.5p . which has been applied throughout and will continue so. We have the RGE analysis 
of this work. We will require both r]i (A) and r^oi (A) to be smaller than vivr, and A (A) to 
be positive. 

There is one important issue to address though before we proceed, and that is how far 
we want the model to be perturbatively predictive and stable. The maximum value Am for 
the mass scale A should not be very high for two reasons. One, more conceptual, is that we 
want to recognize and allow the model to be intermediary between the Standard Model and 
some possible higher structure. The second reason, more practical, is that a too-high A^ 
is too restrictive for the parameters themselves. For example, for the parameters we used 
in the previous sections, in particular mo = 55GeV and mi = 110 GeV, we have seen that 
A gets negative already for A c^ 15TeV whereas rji leaves the perturbativity region much 
later, for A ~ IGOOTeV. The situation can be reversed. For example, for mo = 67GeV, 
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FIG. 8: The running of the Higgs self-couphng A (fuh RGE). It gets negative at about 15TeV for 
this set of parameters' values. 

nil = 135GeV, and 6 = 15°, A can live positive until about 400TeV whereas rji leaves 
perturbativity at about 50TeV. In this section, we set A^ — 40TeV. 

As a third viability tool, we want the model to comply with the measured direct-detection 
upper bounds. In our model, the total cross section for non-relativistic elastic scattering of 



a dark matter WIMP off a nucleon target is given by the relation iJ] : 



C"det 



4:71 {mi\f + mo) f ^ 



,(3) 



(3), 



Xq cos 9 r]Qi sin 9 



-\ 2 



mt 



m^ 



(5.1) 



In this relation, ttijv is the nucleon mass and mB the baryon mass in the chiral limit. The 
quantities Ag and tjq^ are coupling constants of cubic terms in the theory after spontaneous 
breaking of the two symmetries IJ] : 



A 



(3) 



Ao cos 9 + r/oif 1 sin 9] 



(3) 



r/oif 1 cos 9 — Xqv sin 9. 



(5.2) 



The condition we impose is that adet be within the XENON 100 upper-bounds [11 
In work 



15| , we studied phenomenological implications of the model and constraints on 



it, using rare meson decays and Higgs production. A number of inferences were deduced, but 
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FIG. 9: Running of the mutual couplings (full RGE). The inclusion of the other SM particles 
flattens the runnings. 

we will prefer to retain only two. One is that the mixing angle 6 is to be chosen small. This 
is emphasized in view of the mounting evidence of a SM Higgs particle found by ATLAS 
and CMS at the LHC [l|, |2| . The other is that the physical self-coupling Aq is to be small 
too. This was already observed in [l4|, where the relic-density constraint has the tendency 
of 'shutting down' high values of Aq . At the end of the next section, we will comment on 
possible larger values for Aq . 

In this section, the display range of rriQ and mi is from IGeV to 160GeV. Indeed, there 
is no reliable data to discuss regarding a dark-matter mass below the GeV, and in view of 
the behavior of rji at Aq as a function of mi, taking this latter beyond 160GeV is outside the 
perturbativity region. In practice, mo was taken up to 200GeV, with no additional features 
to report. 

Let us start with Aq and 9 both very small. Fig. [12] displays the regions (blue) for which 



the model is viable up to A^ — 40TeV. Here A 



(4) 



0.01 and 6 = 1°. We see that the 



mass mi is confined to the interval llGGeV — 138GeV. The dark-matter mass is confined 
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FIG. 10: Running of the self-couplings (full RGE) with rriQ larger, r/i is more closely tailed by % 
and A decreases and turns negative at about lOTeV. 

mainly to the region above llSGeV, the left boundary of which having a positive slope as mi 
increases. The DM mass mo has also a small showing in the narrow interval 57GeV — 68GeV. 

The effect of increasing the mixing angle 6 is to enrich the existing regions without 
relocating them. This is displayed in Figs. [13] and [14] for which 6 is increased to 5° and 15° 
respectively. We see that, as 6 increases, the region between the narrow band and the larger 
one to the right gets populated. This means more dark-matter masses above 60GeV are 
allowed, but rrii stays in the same interval, roughly 116GeV — 138GeV. 

Increasing the Higgs-DM mutual coupling Aq has the opposite effect, that of shrink- 
ing existing viability regions. Indeed, compare Fig. [TS] for which Aq = 0.1 and 9 = 15° 
(Am — 40TeV) with Fig. [H] We see shrunk regions, pushed downward by a few GeVs, 
which is not a substantial relocation. Remember that increasing Ag raises ?7oi (Aq) well 
enough above 1 so that this latter will leave the perturbativity region sooner. Increasing 
it is also caught up by the relic-density constraint, which tends to shut down such larger 
values of Ag when the dark-matter mass mo is large. The direct-detection constraint has 
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Ao'"^' = 0.01, e = 10°, niQ = lOOGeV, m, = 1 lOGeV 
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FIG. 11: Running of the mutual couplings (full RGE) with ttiq larger, ryoi starts well above Aq and 
Ai and leaves the perturbativity region before the self-coupling rji. 

also a similar effect. 

VI. CONCLUDING REMARKS 

In this work, we have studied the effects and consequences of the renormalization group 
equations at one-loop order on a two-singlet model of cold dark matter that consists in 
extending the Standard Model with two real gauge-singlet scalar fields. The two issues 
we monitored are perturbativity and vacuum stability. The former is controlled by the 
auxiliary-field self-coupling rji and the mutual coupling tjqi between the dark matter and the 
auxiliary fields. The latter is controlled by the Higgs self-coupling A. When the non-Higgs 
SM coupling constants are switched off, all scalar couplings are positive increasing functions 
of the scale A. Reintroducing them flattens the rates for all the scalar couplings and makes 
the Higgs coupling A turn negative at some scale. The mutual couplings Aq (DM-Higgs) and 
Ai (Higgs-auxiliary) stay always well below one, whereas r^oi boosts up for larger ruQ (DM 
mass) and/or mi (auxiliary-field mass), dominating over rji in some regions. 
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FIG. 12: Regions of viability of the two-singlet model (in blue). Physical Higgs self-coupling Ag 
and mixing angle 9 very small. 

We then have investigated the regions in the space of parameters in which the model is vi- 
able. We have plotted these regions in the (mo, mi)-plane while varying the physical mutual 
coupling Aq between the dark matter Sq and the physical Higgs h, and the mixing angle 
6 between h and the physical auxiliary field. We have required that the model reproduces 
the DM relic density abundance, and that it complies with the measured direct-detection 
upper bounds - those of the XENON 100 experiment. We have also imposed the RGE 
perturbativity and vacuum-stability criteria that we deduced from this work together with 
a maximum cutoff A^ — 40TeV, a scale at which heavy degrees of freedom may start to be 
relevant, something that could be probed by future colliders. This analysis has shown that 
for small Aq and 9, the auxiliary-field mass mi is confined to the interval 116GeV — 138GeV, 
while the DM mass rriQ is confined mainly to the region above llSGeV, with a small showing 
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FIG. 13: Regions of viability (blue) of the model. \q still very small, but 9 larger. The region is 
richer, but not relocated. 



in the narrow interval 57GeV — 68GeV. Increasing 6 enriches the existing viability regions 
without relocating them, while increasing Aq has the opposite effect, that of shrinking them 
without substantial relocation. 

It is pertinent at this stage to comment on the implications of the Higgs discovery at the 
LHC on the possibility of having a light dark matter WIMP 5*0 with a mass rriQ < 62GeV, 
a situation allowed in this two-singlet model. Indeed, on the one hand, for such a light dark 
matter, the decay channel h — > SqSq becomes open, and therefore will lower the number of 
Higgs decays into SM particles. On the other hand. The ATLAS and CMS published data 
on Higgs boson searches seem to indicate that the observed boson is SM-like, and so, one 
expects to have stringent constraints on the parameter space when it comes to light dark- 
matter masses. In [l9|, a global fit to the Higgs boson data that includes those presented 
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FIG. 14: The region of viability (blue) is even richer for larger mixing angle 9. 



at the Moriond 2013 conference by the ATLAS and CMS collaborations 20|, l2l|] has been 
performed; see [22] for earlier fits. It has been found that any extra invisible Higgs boson 
decay must be bounded by the following condition on the corresponding branching ratio: 



Br(/i -^ invisible) < 19%. 



(6.1) 



It turns out that in our two-singlet model, the branching fraction of the invisible width of 
the Higgs boson is smaller than the bound above for tuq < 62GeV. Indeed, if we take for 
example niQ = 55GeV used frequently in this work, the ratio T{h -^ SoSo)/T{h — > bb) is less 
than 17%, quite consistent with the above current bound. Therefore, we conclude that the 
two-singlet model is consistent with the current available data regarding the Higgs boson 
searches. 

Finally, we ask whether the model allows for very light cold dark matter. Below 5GeV, 
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An™ = 0.1, e= 15° 



50 



100 



160 



160 



100 



> 

o 



50 



/< 




160 



100 



50 



50 



100 



160 



mo(GeV) 



FIG. 15: The physical Higgs self-couphng Aq shrinks the viabihty region (blue) as it increases. 

direct detection puts no experimental bound on the total cross section cxdet for non-relativistic 
elastic scattering of a dark matter WIMP off a nucleon target. Such a situation allows for 
very small rriQ regions of viability, but only when Aq is quite large (~ 2 and above) and 
not too small ( ~ 15° and above). However, for such values of the parameters, the branching 
fraction of the invisible Higgs decay is larger than 25%, which is excluded by the current 
LHC available data. 
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